High-contrast dispersive readout of a superconducting flux qubit 

using a nonlinear resonator 
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We demonstrate high-contrast state detection of a superconducting flux qubit. Detection is real- 
ized by probing the microwave transmission of a nonlinear resonator, based on a SQUID. Depending 
on the driving strength of the resonator, the detector can be operated in the monostable or the 
bistable mode. The bistable operation combines high-sensitivity with intrinsic latching. The mea- 
sured contrast of Rabi oscillations is as high as 87 %; of the missing 13 %, only 3 % is unaccounted 
for. Experiments involving two consecutive detection pulses are consistent with preparation of the 
qubit state by the first measurement. 

PACS numbers: 03.67.Lx, 85.25.Cp , 85.25.Dq 



Superconducting qubits have been established as 
promising candidates for the implementation of a quan- 
tum information processor . Remarkable achievements 
in this field include the realization of complex single-qubit 
manipulation schemes Q and the generation of entangled 
two qubit states 0,O|. 

For qubit readout, several detectors have been inves- 
tigated experimentally. In general, their efficiency is rel- 
atively poor, for reasons which are presently not well 
understood. Most detectors to date rely on irreversible 
processes in mesoscopic Josephson circuits |fl Bl Isl IflHIoj . 
For these schemes, energy is dissipated on the chip where 
the qubit is placed. Therefore long waiting times are nec- 
essary to bring the qubit, readout, and control circuits to 
their proper initial state. In particular the qubit state is 
strongly disturbed. More recently, dispersive measure- 
ment schemes are bein g in vestigated, which overcome 
these drawbacks [Tl Imlill Hi fl5| . They are based on 
the measurement of the impedance of a resonator coupled 
to the qubit. The energy used to probe this resonator 
is mostly dissipated at a place remote from the qubit 
chip. Also, the backaction on the qubit is low. However, 
the qubit relaxation times are typically comparable to 
the time necessary to for a reliable measurement of the 
impedance. This limits the detection efficiency when a 
linear resonator is used 0, 0] . With a nonlinear res- 
onator, this limitation can be removed, by using a bifur- 
cation transition [lfjj . In this case, the result of a qubit 
measurement is either of two possible oscillation states 
of the driven resonator. These oscillation states can be 
latched and reliably discriminated irrespective of subse- 
quent relaxation of the qubit. A dispersive method using 
latching was successfully demonstrated in [J^, where it 
was used for the readout of a quantronium. 

In this letter we present experimental results on the 
readout of a superconducting flux qubit using a disper- 
sive method. We observed coherent oscillations in two 
distinct operation modes of the detector: monostable for 
weak driving and bistable for strong driving of our non- 
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FIG. 1: Energy (a) and current (b) of the two qubit states 
versus applied magnetic flux, (c) Schematic diagram of the 
qubit and its control and readout circuits. Crosses represent 
Josephson junctions. The probe wave is sent through a bias 
resistor Ri n = 4.7 kfl. The values of the other circuit com- 
ponents are L 3 ~ 0.2 nH, C = 67 pF, M contro i « 15 fH, and 
M readout = 15 pH. 



linear resonator. In the bistable regime the measurement 
contrast is very large, 87 %, which is a significant im- 
provement over previous measurements 0,0 IE EJ EiJ 
. We also performed consecutive measurements of 
the qubit. We find that the results are consistent with 
projection of the qubit state by a first measurement, with 
a probability as high as 90 %. 

In our experiments we use a persistent current qubit 
(PCQ) ^?|- It is formed of a superconducting ring in- 
terrupted by three Josephson junctions. For suitable pa- 
rameters of the Josephson junctions and for an applied 
magnetic flux in the qubit loop $ q b close to (n + 1 /2)<&o 
(n is integer and $0 — h/2e), this system can be operated 
as a two level system. The qubit Hamiltonian is 



H qb = \(eu z + Aa x 



(1) 



where e = 2/ p $o/- A and I p are fixed parameters, de- 
termined by the charging and Josephson energies of the 
Josephson junctions, and the frustration / = (& q b — (n + 
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1 /2)4> )/ <J?o- Figure^, shows a plot of the ground (g) 
and excited (e) energy levels of the qubit used in our 
measurements. The expectation value of the circulating 
current in the qubit ring corresponding to each eigenstate 
is plotted in Fig. ^p. 

Figure^ shows a schematic representation of our read- 
out and control circuit. The readout circuit is based on 
a DC-SQUID, coupled to the qubit through their mu- 
tual inductance M rea( i out . The SQUID has a Joseph- 
son inductance 0] Lj($ sq ), which depends on the total 
magnetic flux in its loop $ sq = $ sqM as + M readout I qb ; 
&sq,bias and M rea d utlqb are generated by an external coil 
and by the qubit, respectively. The contribution of the 
self-generated flux is negligible compared to & S q,bias and 
nearly constant over the range of variation of & S q,bias 
relevant in our experiments. The two qubit states have 
different currents I q b (see Fig. ^>) , resulting in values of 
the Josephson inductance different by « 2 % in our ex- 
periment. We form a resonant circuit by including Lj 
in a loop closed by the capacitors C and the small un- 
intended stray inductance L s (see Fig. We probe 
the transmission of this resonant circuit by sending a 
wave with a frequency F and current amplitude at the 
sample I r ,o and measuring the amplitude V r ^o and rel- 
ative phase x °f the voltage of the transmitted wave. 
We use a low- noise cryogenic amplifier, described in [l9j . 
When F is close to the resonance frequency F res , the 
transmission depends strongly on small changes of Lj 
and thus on the qubit state. In our experiment F res ~ 
800 MHz and the quality factor Q w 100. The qubit 
state is controlled by applying a magnetic flux in its ring 
® q b(t) = <b q b,bias + M contro iI contro i(t). $ q b,bias is gener- 
ated using an external coil and set to a value close to one 
of the operation points (n + l/2)$ - M contro il control (t) 
is generated by an on-chip line and used to change the 
qubit state either adiabatically or by resonant pulses. 
The qubit, resonant circuit, and control line are fabri- 
cated on an oxidized silicon substrate by using standard 
electron beam lithography. Measurements are performed 
in a dilution refrigerator, at temperatures of rs 30 mK. 

The Josephson inductance of the SQUID depends on 
the current in a nonlinear way. Our resonant circuit is de- 
scribed, to a good level of approximation, by the Duffing 
oscillator model |2jj. This model predicts that a certain 
driving condition of the resonator (given by the driving 
frequency F and amplitude I r o) results in one unique 
forced oscillation state when F > F res (l — \/3/(2Q)). 
When F < F res (l — >/3/(2Q)), two oscillations states 
are possible for Ibi < Fa < Ibh, and only one state 
otherwise. Ibi and Ibh are the values of the driving cur- 
rent where bifurcation occurs, dependent on F and the 
parameters of the resonant circuit. The oscillation state 
denoted by exists when I t q < Ibh-, and the state 1 when 
I r o > Ibi- Without fluctuations the resonator, driven at 
frequency F and amplitude I r o slowly increasing from 
IrO =0; resides in state until I r o reaches Ibh-, when it 
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FIG. 2: Measurements in the monostable regime. (a) 
Schematic representation of the qubit control (top) and read- 
out (bottom) sequence; explanations are given in the text, 
(b) Spectroscopy measurement for f m = 0, f p = f r = 0.0021, 
T m = 2.4 us, and T r = 6 /is. (c) Plot of the qubit transition 
frequency, measured spectroscopically, versus the frustration 
f m - The solid line is a fit that yields the qubit parameters, 
A and l p , given in the text, (d) Measurement of Rabi oscil- 
lations with the same parameters as in (c) and F m — 4.602 
GHz, resonant to the qubit energy level splitting. The solid 
line is a fit with an exponentially damped sinusoidal. 



switches to state 1. Thermal fluctuations will cause the 
transition — > 1 to occur randomly at J r o < Ibh, with a 
rate that increases exponentially as I r o approaches Ibk- 
The bistable measurement regime relies on the fact that 
the bifurcation current Ibh depends on the value of the 
magnetic flux in the SQUID loop $> sq . For a fixed driving 
condition of the oscillator, the probability for the transi- 
tion — > 1 to take place is very sensitive to variations in 
& sq , and thus to changes in the qubit state. The monos- 
table measurement regime does not exploit the — ► I 
transition, but relies instead on the measurement of the 
voltage and phase of the driven resonator, which also de- 
pend on Q sq . 

We present our experimental results starting with mea- 
surements in the monostable regime. We control the 
qubit state by using a pulse as shown in Fig. |2K (top). 
We start with a preparation part (P), by fixing the frus- 
tration / to a value f p , for a time much longer than the 
qubit energy relaxation time. This results in a Boltz- 
mann distribution for the qubit initial state, set by the 
effective temperature T qo ^ e ff. The qubit is prepared in 
the ground state if T q b } eff <C E e — E g . After qubit state 
preparation, the frustration / is changed adiabatically to 
the value f m , where qubit manipulation (M) is done by 
using a microwave burst of frequency F m and duration 
T m . Finally, the frustration / is changed adiabatically to 
the value f r (for readout - R), preserving the weight of 
the energy eigenstates created at M. Readout is possible, 
provided e(/ r ) > A, because the two qubit states have 
significantly different currents (see Fig. \T]p) . To measure 
the state of the qubit, a microwave burst is sent to the 
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resonator with frequency F, amplitude I r o,max, and du- 
ration T r (see Fig. [SJi-bottom). The sequence shown in 
Fig. [2t is repeated typically 10 5 times, and the average 
of the voltage of the transmitted wave (V r o), which is 
a measure of qubit population, is calculated. Figure [2b 
shows the result of a typical spectroscopy measurement. 
The position of the qubit spectroscopy peaks, F pea k, is 
plotted as a function of the frustration f m in Fig. [2t. A 
fit of the data with the expectation for the qubit energy 
level splitting following from Eq. ^ yields the qubit pa- 
rameters A = 4.58 GHz and I p = 360 nA. Using strong 
microwave pulses, with a frequency F m equal to the en- 
ergy levels splitting of the qubit, we induce Rabi oscil- 
lations between the qubit energy eigenstates. A typical 
result is shown in Fig. Eli, together with a fit with a ex- 
ponentially damped sinusoidal signal. We observe Rabi 
oscillations with frequencies as high as 100 MHz, varying 
linearly with the amplitude Sf m of the microwave burst. 

We now turn to a presentation of the experiments us- 
ing the bistable operation mode. To characterize experi- 
mentally the transition between the states and 1 of the 
driven oscillator, we use the methods developed in [l6| . 
We drive the resonator with an AC pulse of fixed fre- 
quency F < F res (l — v3/(2Q)) and modulation of the 
amplitude as shown in Fig. [3]}. During the first part of 
the pulse, of duration T SWl the amplitude is set to a value 
Iro.sw close to Ish, such that there is a significant proba- 
bility for the — > 1 transition of the oscillator. The sec- 
ond part of the pulse is used to discriminate the states 
and 1 by measuring the amplitude and phase of the 
voltage of the resonator. In this part I r o is set such that 
both transitions — ► 1 and 1 — > have negligible rates 
and a long time Ti atc h allows for amplifier noise suppres- 
sion. We first measure the probability P sw of switching 
between the states and 1, with the qubit in its ground 
state, for two values of the magnetic flux in the SQUID 
loop QsqMasi and & S q,bias2 ■ Their difference equals the 
change in magnetic flux when the qubit state changes 
from g to e. The switching probability P sw is shown in 
Fig- EI 3 as a function of the amplitude I r o,sw The figure 
shows that it is possible, in principle, to detect the qubit 
state with an efficiency as high as 98 %. 

Using the bistable operation mode, we next measure 
changes in the state of the qubit. We use a control pulse 
as shown in Fig. |2t (top) and a readout pulse as shown 
in Fig. The shape of the control pulse is optimized 
to achieve high fidelity for qubit state preparation. Fig- 
ure (SJ; shows the switching probability curves when the 
qubit is prepared in the g or e state by using a resonant 
microwave pulse at M (see Fig.|2K)- The maximum sep- 
aration between the two curves is 87 %. This is smaller 
than the value of 98 % corresponding to a measurement 
of a magnetic flux difference equal to the qubit signal 
(see Fig. Ob). Detailed experimental studies indicate that 
the most important contribution to the observed loss of 
contrast is the adiabatic shift from M to R, accounting 
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FIG. 3: Measurements in the bistable regime, (a) Envelope 
of the AC readout pulse, (b) Switching probability curves for 
two values of $ s ,,w„ s : 4> S9 ,bm S i = 2.335 $o and <5> sqi u aS 2 = 
2.340 $o- (c) Switching probability curves, for qubit in the g 
and e states, for f p = 0.0161, f m = 0.00, f r = 0.0126, T sw » 
100 ns, and Tiatch ~ 10 /is. (d) Rabi oscillation with the same 
setting as for (c), and I 3W optimized for maximum contrast. 



for « 7 %. The qubit relaxation during the measure- 
ment (but not due to the measurement process) accounts 
for w 2 %. The remainder of « 3 % occurs during the 
switching on of the detector. In Fig. [St we show the mea- 
surement of Rabi oscillations with measurement settings 
optimized for maximum contrast. Each measurement is 
repeated ~ 10 4 times in order to reduce the spread of the 
average when measuring superpositions. 

We finally address the question of how the qubit state 
is changed by the measurement |2lj . This requires the 
use of two consecutive measurements; the second mea- 
surement is used to characterize the effect that the first 
measurement had on the qubit. Using two consecutive 
readout pulses, of the type shown in Fig. has the 
disadvantage that the minimum duration of the latch- 
ing plateau, due to the noise of our amplifiers, is com- 
parable to the qubit relaxation time. The qubit would 
undergo significant energy relaxation before the second 
measurement takes place. We choose instead a readout 
pulse as shown in Fig. 0k, consisting of two intervals 
of large amplitude, realizing two consecutive measure- 
ments Ri and R2. We use the labels n and r 2 to de- 
note the outcome of the measurements Ri and R2; they 
take the values (1) if switching did not (did) occur. 
For this double-measurement pulse the probability that 
switching occurs is P^ w (\ip) in it, Ir0, sw > I? , sw , T d eia V ) 
I — P(n = A r 2 = 0); here \ip)i n u is the initial 
qubit state, J* sw and Iro,sw are t ne amplitudes for Ri 
and R2, 1 'delay is the delay between the two measure- 
ments, and P(r\ = A r 2 = 0) is the probability that 
switching does not occur during either Ri or R 2 . We 
also read out the qubit state with a single measurement 
pulse by setting the amplitude I^ sw = I r o,iatch- 

The 

switching probability for this single-measurement pulse is 
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FIG. 4: Results of the correlation measurements, (a) En- 
velope of the double-measurement pulse, (b) The switching 
probability curves for the qubit in the ground or excites states, 
with a single-measurement pulse with T$ m m 125 ns. (c) The 
contrast of Rabi oscillations measured with a pulse as shown 
in (a), with I^ S w = =7.7 nA and T} w = T% w » 125 ns, 

as a function of Tdeiay The solid line is a fit with an expo- 
nential decay, (d) Conditional probability for the qubit to be 
in the excited state, when the first measurement did not lead 
to switching, for different qubit initial states; Tdeiay = 800 ns. 



PLMinit,% , sw ) = 1 - P(n = 0), where P(n = 0) is 
the probability that switching does not occur during the 
measurement Ri. Figure 0Jd shows the switching proba- 
bility curves, for a single-measurement pulse, as a func- 
tion of Iro sw - (The timing for this single-measurement 
pulse is as used in the experiments described below, but 
different from those that led to the data shown in Fig.^J:). 

In a first experiment we measure the contrast of Rabi 
oscillations with a double-measurement pulse as a func- 
tion of the delay time Tdeiay between Ri and R2, as shown 
in Fig.^Ji. The observed decrease is due to the relaxation 
of the qubit between Ri and R2 . A fit with an exponen- 
tial decay, shown by the solid line in Fig.QJ:, yields a de- 
cay time which is the same as the qubit relaxation time. 
This shows that the correlations between the results of 
measurements Ri and R2 are due to the change in the 
qubit state, and not to the detector dynamics. 

In a second experiment we set the initial state of the 
qubit to a defined value and we measure the joint proba- 
bility P(n = A T2 = 0) with a double-measurement 
pulse and the probability P(r± = 0) with a single- 
measurement pulse. This allows calculating the prob- 
ability for T2 — when r\ = 0. Therefore the con- 
ditional probability P(q — e\r\ — 0), for the qubit 
to be in state e when the first measurement gave the 
result r\ — 0, can be calculated. This is plotted in 
Fig. Hi versus i^Os-ici f° r initial qubit states given by 
W)imt — cos(0/2)|<7) + exp(i(f>) sin(0/2)|e) obtained by 
Rabi rotation; 8 is indicated in Fig. 0Ji for each curve, 
and 4> is an unimportant phase angle. For large sw 
all the curves collapse to a single value of about w 90 %. 



A no-switch event results in a post-measurement state 
nearly equal to the excited state, irrespective of the pre- 
measurement state. The region over which the curves 
in Fig. 0Ji show a fast increase corresponds to the re- 
gion where the switching probability for the ground state, 
shown in Fig. 0Jd, increases; this can be interpreted as a 
transition from weak to strong measurement. Note that 
using the double- measurement pulse shown in Fig. Qli in- 
stead of two successive pulses as shown in Fig.|2K has the 
limitation of only giving information on the qubit state 
when n = 0. 

In this paper we presented measurements of a super- 
conducting qubit using a dispersive readout scheme. Op- 
eration in the bistable mode results in the observation 
of Rabi oscillations with very high contrast, significantly 
improved over other measurements of this type. We also 
present measurements which indicate that a single mea- 
surement prepares the qubit state with a fidelity of about 
90 %. These results establish dispersive readout as very 
suited for flux qubits and make it promising for readout 
of entangled two-qubit states. 
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